CLOSED TRAJECTORIES OF A PARTICLE MODEL 
ON NULL CURVES IN ANTI-DE SITTER 3-SPACE 



EMILIO MUSSO AND LORENZO NICOLODI 

Abstract. We study the existence of closed trajectories of a particle 
model on null curves in anti-de Sitter 3-space denned by a functional 
which is linear in the curvature of the particle path. Explicit expressions 
for the trajectories are found and the existence of infinitely many closed 
trajectories is proved. 



1. Introduction 

In this paper we study null curves in anti-de Sitter 3-space which are 
critical points for the functional 

(1) £(7)= / (m + hy)ds, mel, 

J-y 

where s is the pseudo-arc parameter which normalizes the derivative of the 
tangent vector field of 7 and £7(5) is a curvature function that, in general, 
uniquely determines 7 up to Lorentz transformations ([1], [S]). The func- 
tional ([!]) is invariant under the group SL(2,R) x SL(2,M), which doubly 
covers the identity component of the group of Lorentz transformations. Mo- 
tivations for this study are provided by optimal control theory and especially 
by the recent interest in certain particle models on null curves in Lorentzian 
3-space forms associated with action integrals of the type above ([11]. [1]). 
Yet another motivation is given by surface geometry; if we take SL(2,R) as 
a model for anti-de Sitter 3-space, then a null curve in anti-de Sitter 3-space, 
as real form of a holomorphic null curve in SL(2, C), is related to the theory 
of constant mean curvature one (cmc-1) surfaces and flat fronts in hyper- 
bolic 3-space ([2], [13], [5]). In perspective, one would like to understand 
the class of cmc-1 surfaces generated by the critical points of (pQ). 

The purpose of this article is to investigate the global behavior of extremal 
trajectories of the functional ([I]). We will find explicit expressions for the 
extremal trajectories and then establish the existence of infinitely many 
closed ones. This result is related to the presence of maximal compact 
abelian subgroups in the isometry group SL(2, M) x SL(2, R) of anti-de Sitter 
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3-space. For a discussion of extremal trajectories in the other Lorentzian 
space forms we refer to [6], [10] . 

The Euler-Lagrange equation associated to ([I]) yields that the curvature 
of extremal trajectories is either a constant, in which case we have null 
helices, or an elliptic function^ of the pseudo-arc parameter. In this case, 
extremal trajectories are governed by a second order linear ODE with doubly 
periodic coefficients. By classical results of Picard [12] in the Fuchsian theory 
of linear ODEs, the solution curves are then expressible in terms of the 
Weierstrass p, a and £ functions (cf. Theorem [1]). The explicit integration 
of extremal trajectories amounts to the integration of a linearizable flow on 
K\ x K2, where K\ and K2 are 1-dimensional abelian subgroups of SL(2,R). 
In particular, if K\ = Ki = SO(2), the integration amounts to solving a 
linearizable first-order ODE on a 2-dimensional torus. This setting strongly 
suggests the possibility of periodic solutions. That this is indeed the case 
is established in Theorem [2] where the existence of countably many periodic 
trajectories is proved by studying the map of periods. The proof relies on 
computations made with Mathematica. 

The paper is organized as follows. Section [2] contains some background 
material. Section [3] provides the explicit integration of extremal trajectories. 
Section [4] discusses periodic trajectories and proves the existence of infinitely 
many of them. The Appendix outlines the derivation of the Euler-Lagrange 
equation associated with (pQ) via the Griffiths formalism [7j. 



Anti-de Sitter 3-space, Wf, can be viewed as the special linear group 
SL(2, R) endowed with the bi-invariant Lorentz metric of constant sectional 
curvature —1 defined by the quadratic form 



for each X = {x)) £ s((2,R). The group G = SL(2,R) x SL(2,R) acts 
transitively by isometries on Hf = SL(2,R) via the action 



2. Preliminaries 



(2) 



q{X) = {x\f + x\x\ 



det X, 




tt: G 3 (A,B) i-» AB- 1 e SL(2,R) 



makes G into a principal bundle with structure group A. 



Possibly a degenerate one, i.e., an hyperbolic, trigonometric, or rational function. 
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Let JcRbe any open interval of real numbers. A smooth parametrized 
curve 7 : I — > SL(2,R) is null, or lightlike, if det (7 _1 7') vanishes identically. 
If 7 has no flex points0 there exists a canonical lift 

r = (r+,r_) : j-g 

such that 

where a; is a nowhere vanishing 1-form, the canonical arc element, and A; : 
I — ► R is a smooth function, the curvature function. We call T the spinor 
frame field along 7 and its components T + and T_ the positive and negative 
spinor frame, respectively. The spinor frame T is essentially unique, in the 
sense that ±r are the only lifts satisfying ([3]). Throughout the paper we 
will consider null curves without flex points and parametrized by the natural 
parameter, i.e., to = ds (cf. [3J, [3]). 

Conversely, for a smooth function k : I — ► R, let H±(k) : / — ► st(2,R) be 

<") H +w = ( ft +i J)' "- w = (<=°i J 

By solving a linear system of ODEs, there exists a unique (up to left multi- 
plication) 

r = (r+,r_) -.i^g 

such that 

T- 1 T' + = H + (k), TZ 1 T'_ = H_(k). 

In particular, 7 = T+n 1 : I — > SL(2,R) is a null curve without flex points, 
parametrized by the natural parameter and with curvature function 

In this context, two null curves 7 : J — > SL(2,R), 7 : J — > SL(2,R) 
are said equivalent if there exist c € R and ^,5 6 SL(2,R) such that 
7(s) = A^(s + c)S _1 , for all s G /. 

Definition. An extremal trajectory (or simply, a trajectory) in is a 
null curve with non-constat curvature which is a critical point of the action 
functional 

(5) C m (l) = / (m + k^u 

under compactly supported variations, where the Lagrange multiplier m is 
a real constant. 

The Euler-Lagrange equation associated to C m is computed to be 

k'" - 6kk' + 2mk' = 



2 i.e., y'(t) and y"(t) are linearly independent, for each t € I, where 7" is the covariant 
derivative of 7' along the curve. 

The curve 7 is uniquely defined up to orientation and time-orientation preserving 
isometries. 
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(cf. [J| ; see also the Appendix for a different way of deriving this equation) . 
This may be thought of as the intrinsic equation of a trajectory. If we let 
h := 5 (A; — y) be the reduced curvature, then h satisfies 

(6) (h'f = 4h 3 -g 2 h - 53) 

for real constants g 2 and g 3 . Hence h is expressed by the real values of either 
a Weierstrass p-function with invariants g 2 , 93, or one of its degenerate 
forms. 

We call a solution to ([6]) a potential with analytic invariants g 2 , 53- Two 
potentials are considered equivalent if they differ by a reparametrization 
of the form s •— ► s + c, where c is a constant For real 52 and 53, let 
A (32 j 03) = 27g 3 — g\ be the discriminant of the cubic polynomial 

P{t;g 2 ,g 3 ) = 4t 3 -g 2 t-g 3 - 

The study of the real values of the Weierstrass p-function with real invariants 
92 , 93 (and its degenerate forms) leads to primitive half-periods ui\ , uj 3 such 
that: 

• A(g 2 ,g 3 ) < 0: > 0, uj 3 = ivui, v > 0. 

• A(g 2 ,g 3 ) > 0: u\ > 0, uj 3 = \{l + iv)ui, v > 0. 

• ^(92,93) = and g 3 > 0: Ui > 0, u; 3 = -Moo. 

• A(g 2 ,g 3 ) = and 53 < 0: u)\ = +00, -iu> 3 > 0. 

• 92 = 93 = 0: Wi = +00, uj 3 = +ioo. 

Accordingly, denoting by T>(g 2 ,g 3 ) the fundamental period-parallelogram 
spanned by 2u\ and 2co 3 , the only possible cases for the potential function 
h : I -> R are: 



/i(s) = p(s;5 2 ,93), ^ = (0,2o;i). 
M s ) = P3{s;g2,93) = p(s + U3;92,93), I 
h(s) = p(s;g 2 ,g 3 ), I = (0,2ui). 



A < 
A < 
A > 

A = 0, g 3 = -8 a 3 > 



h(s) = —3a tan 2 (v 7 — 3as) — 2a, I = (- 



A = 0, 53 = -8a 3 < 0: 

a(s) = 3atanh 2 (v 7 ^) - 2a, I = R. 
52 = 93 = 0: = s~ 2 , I = (-00, 0) or / = (0, +cxd). 



4 When invariants (?2 and gs are given, such that 27<?f ^ pi, the general solution of the 
differential equation (^f ) 2 = 4y 3 — giy — <?3 can be written in the form p(z + a; (72,53), 
where a is a constant of integration. 
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3. Integration of the extremal trajectories 
For a potential function h with invariants <?2, 93, let 



H±{m,h) := \\j p (y ± !; 52,53 ) 
and, for each so £ define (j)±(m, h, sq) : I — ► R by 



4>±(m,h,s ) 



h{u) 



1 



±1) 



d«, n±(m, h) 7^ 0, 
/J>±(m, h) = 0. 



fc(u) - (f ± 1) 

Next, let ui±(m, /i) be the unique point^l in V{g2,gz) such that 



h(w±) = — ±l and h'(w±) = 2fj,±(m, h). 

Then, denoting by an and 0i> respectively, the sigma and zeta Weier- 
strassian functions corresponding to the potential /ijl we compute: 

Case I: if /x±(m, h) / 0, 

</>±(m, h, s ) = log ah [ S W± \ + 2s( h (w±) + c(s ). 



Case II: if fi±(m, h) = and p| + ffl / 0, 

-1 



(f)±(m,h,s ) 



3(f ±l) 2 - 52 /4 



{c*( 



so)- 



Case III: if //±(m, /i) = and 52 = 53 = 0, 

(f)±(m,h,s ) = ^s 3 + c(s ). 

Accordingly, define the maps R±(m,h), D±(m, h, so) : / — ► GL(2, C) as 
follows: 



Case I: if n±(m, h) ^ 0, 



/ /i' — 2fi±(m, h) 



R±(m.h) 



2fi± 



f ±1) 



D±(m,h,s ) 



2Jh 

h' + 2n±(m,h) 
exp(-0 ± (m,/i,s o )) 



^-(f ±1) 



h-(?±l) 





exp ((^>±(m, /i, so)) 



If m = ±3 and 32 = 33 = 0, u>± = 00. 

0^ and £h are the unique analytic odd functions whose meromorphic extensions satisfy 
Ch = -ft and a' h /a h = < h 
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Case II: if /i±(m, h) = 
R±(m.h) 



( 1 

h! 



[f ±1) 



(j)±(m,h,so) 



D±(m,h,a ) = \ y/h-(f±l 

1 

Finally, define the maps F±(m, h, sq), 7(771, h, sq) : I — ► GL(2, C) by 
r ± (m, h, s ) := R±(m, /i)(s ) _1 .D±(m, h, s )(s )' 1 D±(m, h, s )R±(m, h), 
7(771, /i, s ) := r + (m, /i, s ) r_(m, /i, so) -1 . 
We are now in a position to state the following. 

Theorem 1. The curve 7(771, h, sq) takes values in = SL(2,R) and de- 
fines an extremal trajectory with multiplier m and reduced curvature h. In 
particular, any extremal trajectory is equivalent to a curve of this type. 

Proof. A direct, lengthy computation shows that 

(7) r±(m, h, s ) _1 r±(m, h, s )' = H(m, h)±, 
where 

(8) H + (m,h)=( 2h + ° f + 1 J), JT_(m,fc)=( 2/i+ ° f _ 1 J 

Equations ([7]) and (J8|) imply that T±(m, h, sq) take values in a left coset of 
SL(2,R). Our normalization implies that T±(m, h, sq)(sq) = I2, and hence 
that T±(m, h, sq) and 7(771, h, sq) take values in SL(2,R). On the other 
hand, ([7|) and (|8|) imply that T(m, h, so) = (T + (m,h,so),F-(m,h,so)) is a 
spinor frame field with multiplier m and curvature function k = 2h + y. 
Therefore, 7(771, /"i, so) is a trajectory with curvature function k. This proves 
the required result. □ 

4. Periodic trajectories 

A trajectory is said quasi-periodic if its reduced curvature is a periodic 
function and if ^±{m,h) are purely imaginary. Let Q (t;£, e%) be the cubic 
polynomial 

The reduced curvature of a quasi-periodic trajectory can be written in the 
form 

(9) h{s;£,ex) = e\ 
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where (m,£, e±) belongs to 




) <0} 



The period of h(s; £, e\) is given by 



2 - P 

p(e, ei ) = 2J^—K(£), 



where K is the complete elliptic integral of the third kind. We put 



where C G GL(2, C) and M : M. — > GL(2,C) is a periodic map with period 
p(£, ei). If we set H(m, £, ei) = V(p(£, ei),m, £, ei), for each (m, £, e\) G W, 
then a quasi-periodic trajectory with invariants (m, £, e\) G W is periodic if 
and only if LT(m, £,ei) £QxQ. 
We can state the following. 

Theorem 2. There exists a discrete subset M.q C R suc/i /or every 
m G ]R\A4o, i/iere exisi countably many periodic trajectories with multiplier 
m. 




and let V = (V+,V-) : K x W -> M 2 be the analytic map 




Recall that the spinor frame fields of h are of the form 




Figure 1. The graph of the function /. 



Proof. Consider the analytic map 



* : W 3 (m,£,ei) i-> 



d(n + ,n_) 

0(4 ei) 



l(m^ei) £ R- 
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If W(m,£,ei) + 0, then U m : {i,e x ) -> U(m,£,ei) 6 K 2 is a local dif- 
feomorphism near (£, el). In this case there exist countably many closed 
trajectories with multiplier m. The mapping ^ can be computed explicitly, 
or numerically. To avoid dealing with quite long formulae, we would rather 
adopt the numerical viewpoint. Once we know we define / : K — ► R, 
m i — > A00^(m, i, \m\ + 10). This is an analytic function for m < and for 
m > 0. Looking at the graph of / (cf. Figure 1), we see that / 7^ const. 
This implies that the set Mq = {m E R| f(m) = 0} is a discrete set and 
that, for every m £ R \ .Mo, there exist countably many closed trajectories 
with multiplier m. □ 

5. Appendix: Derivation of the Euler-Lagrange equation 

In this section, we outline the derivation of the Euler-Lagrange equation 
associated with the G- invariant functional (pQ). We follow a general con- 
struction for invariant variational problems with one independent variable 
(due to Griffiths) and write the Euler-Lagrange equation as a Pfaffian dif- 
ferential system (PDF) {J , uj) on an associated manifold Y. We adhere to 
the terminology and notations used in [7j. 

5.1. The variational problem. The starting point of the construction is 
the replacement of the original variational problem on null curves in SL(2, R) 
by a G-invariant variational problem for integral curves of a Pfaffian differ- 
ential system (PDF) with an independence condition on M := G x R. 
Let (a, (3) be the Maurer-Cartan form of G = SL(2, R) x SL(2, R), where 

1 - 1 X 'ft d\ 



(10) a ={ 4 13 ~ \ ft 2 -A 

The Maurer-Cartan equations of G, or the structure equations, are given by 

'dl3\ = -(3l/\aj 



da\ = —a\ A a\ 
da\ = 2a\ A a\ 
da\ = —2a\ A a\ 



df3\ = -20\ha\. 



On M := G x R, consider the PDS (X, uj) defined by the differential ideal 
X generated by the linearly independent 1-forms 



2V«2 



01), V 2 = a{, V 3 = 0l, 



where 



l(a\ + 0l)-u, rf = \{a\+0l)-ku, 

u ■= - 0l) 



gives the independence condition uj 7^ 0. Now, if 7 : I — > SL(2, R) is a 
null curve without flex points, then the curve (r 7 ;A; 7 ) : / — > M, whose 
components are, respectively, the spinor frame field T = (r + ,r_) along 
7 and the curvature of 7, is an integral curve of (T, uj) (cf. Section [2]). 
Conversely, any integral curve (r, k) : I — > M of (T, uj) defines a null curve 
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with no flex points 7 : I — > SL(2, M), where T is the spinor field along 7, and 
k is the curvature of 7. So, null curves without flex points in SL(2,M) are 
identified with the integral curves of (I,oj). 

From the above discussion, it follows that a null curve 7 without flex 
points is an extremal trajectory if and only if the pair (T 7 , /c 7 ) of its spinor 
frame field T 7 and curvature /c 7 is a critical point of the functional defined 
on the space V(T, uj) of integral curves of (I, uj) by 

(11) C : (r, k) G V(2, u) j {T,k)*((m + k)u) 

when one considers compactly supported variations. Here, I(r,k) is the do- 
main of definition of (T, k). 



5.2. The Euler Lagrange system. Following Griffiths [7], the next step 
is to associate to the variational problem (jl 11) a PDS {J ', uj) on a new man- 
ifold Y, whose integral curves are stationary for the associated functional. 
For this, let Z C T*M be the affine subbundle defined by 

Z = (m + k)u + / C T*M, 

where / is the subbundle of T*M associated to the differential ideal I. The 
1-forms (77 , . . . , rf , to) induce a global affine trivialization of Z, which may 
be identified with M x R 5 by setting 

M x M 5 3 ((T,k);xi,. . . ,x 5 ) ^ uj mk) + xjrf l{m G Z 

(we use summation convention). Accordingly, the Liouville (canonical) 1- 
form of T*M restricted to Z is given by 

H = (m + k)uj + Xjrf . 

Exterior differentiation and use of the structure equations give 

dfi = dk A uj + (m + k) [(1 + k)r] 2 + (1 - k)rf] Auj + dxj A rf 
+ xi{rf - rf) A uj + x 2 [rf - (1 + k){rf + r/ 4 )] A uj 
+ ^3 [rf + {k- l){rf - r] 4 )] Au-14 [(2 + k)rf + (2 - k)rf] A uj 
+ x 5 [(1 - k 2 )rf - (1 - k 2 )rf )\ Auj - x 5 dk A a; mod {rf A }. 



Then, we compute the Cartan system C{dii) C T*Z determined by the 2- 
form d[i, i.e., the PDF generated by the 1-forms {i^dii \ £ G %(Z)} C fi^Z). 

Contracting d^i with the vector fields of the tangent frame ^^j, Jr, 

on Z, dual to the coframe (uj,dk,r]\dxj), i,j = 1,...,5, we establish the 
following. 
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Lemma 3. The Cartan system C(d(j), with independence condition uj, is 
generated by the 1- forms r/ , . . . rf and 

7Ti = (1 — x 5 )dk, 7r 2 = (1 — x 5 )lo, 0i = dxi + [x 2 (l + k) + x 3 (l - k)] U, 
fa = dx 2 + [a?i - (m + fc)(l + k) + x±(k + 2) + x 5 (k 2 - 1)] w, 
/% = dx 3 + [-xi — (m + k)(l - k) + x 4 (2 - k) + x 5 (l - A; 2 )] w, 
k /?4 = <ix 4 + [x 2 (l + fe) + x 3 (/c - 1)] u, f3 5 = dx 5 - (x 2 + x 3 )uj. 

The Euler-Lagrange system associated to the variational problem is the 
PDS {J ,uj) on a submanifold Y C Z obtained by computing the involutive 
prolongation of (C(dfj,),u). The submanifold Y is called the momentum 
space. A direct calculation gives the following. 

Lemma 4. The momentum space Y C Z is defined by the equations 

m + k 

£5 = 1, x 2 + x 3 = 0, x 4 = — - — . 

The Euler-Lagrange system (J, to) is the PDF on Y with independence con- 
dition lj generated by the 1- forms n 1 ^, ■ ■ ■ , ?? 5 |y and 

f u\ = dx\ + 2kx2UJ, 

k 2 mk 



cr 2 = dx 2 + y- 1 + xi j lu, 

0"3 = dk + \X20J- 

Remark 5. The importance of this construction is that the projection 7ry : 
Y — ► M maps integral curves of the Euler-Lagrange system to extremals 
of the variational problem associated to (M,X). In our case, the converse 
is also true (see also below), so that all extremals arise as projections of 
integral curves of the Euler-Lagrange system. The theoretical reason for 
this is that all derived systems of J have constant rank (cf. [I], [8]). 

A direct calculation shows that //iy A(d/U|y) 4 7^ on Y, i.e., the variational 
problem is nondegenerateQ This implies that /xiy is a contact form and 
that there exists a unique vector field £ £ 3£(Y"), the characteristic vector 
field of the contact structure, such that Miy(C) = 1 and i^dpLvy = 0. In 
particular, the integral curves of the Euler-Lagrange system coincide with 
the characteristic curves of £. 

5.3. The Euler-Lagrange equation. Let V(J,uj) be the set of integral 
curves of the Euler-Lagrange system. If y = ((T, k); zi,a; 2 ) : I — ► K is in 
V( 1 7,a;), the equations 

n 1 = n 2 = ■ ■ ■ = 7] 5 = 



'''A variational problem is said to be nondegenerate in case 

dimY = 2m + l and /i|y A (d/i|y) m 7^ 0. 
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and the independence condition tell us that T = (r+,r_) defines a 

spinor frame along the null curve 7 = r+ri 1 and that k is the curvature of 

7- 

Next, for the smooth function k : I — ► R, define fe', fe" and k'" by 
dfc = k'u, dk' = k"u, dk" = k"'oo. 
Equation 03 = implies 

y 



X2 = ~J 



Further, equation a 2 = gives 



k" k 2 mk 

Finally, equation a± = yields 

(12) k'" - 6kk' + 2mk' = 0. 

This coincides with the Euler-Lagrange equation of the extremals of ([I]), 
which has been computed in [5]. Thus, an integral curve of the Euler- 
Lagrange system projects to an extremal trajectory in SL(2,R). 

Conversely, if 7 : I — > SL(2, R) is a null curve without flex points, T 7 its 
spinor frame, and k^ its curvature, let y 7 : I — > Y be the lift of 7 to Y given 
by 

k" k 2 mk _ k n 



y 7 (t) = I (r 7 , fey); x~T + ~ + 1 '~T 

Then, y 7 is an integral curve of the Euler-Lagrange system if and only if fc 7 
satisfies equation (|12p if and only if 7 is an extremal trajectory. Thus, the 
integral curves of the Euler-Lagrange system arise as lifts of trajectories in 
SL(2,R). 

Remark 6. Griffiths' approach to calculus of variations, besides for providing 
the Euler-Lagrange equations, is important for giving an effective procedure 
to construct the momentum mapping induced by the Hamiltonian action of 
G on Y and to prove that it is constant on the integral curves of (jT, u), 
which in turn leads to the integration by quadratures of the extrema (cf. 
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